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Abstract 

Motivated by the recent BaBar measurements of B — > D^*H~ui semileptonic decays and the 
resulted 3.4a discrepancy of the measured values of R(D^*') from the standard model (SM) 
predictions, we calculate the branching ratios B(B — > D^*H~i>i) and the ratios R(D^) by 
employing the perturbative QCD (pQCD) factorization approach. We find that the pQCD 
predictions are R{D) = 0.429±0.023 and R(D*) = O.SOl^;^: the BaBar's discrepancy of 2.0a 
for R(D) and 2.7a for R(D*) are therefore reduced significantly to 0.2a and I.Oct respectively. 
We provide successfully the first SM interpretation for the BaBar's R(D^) anomaly. 
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Very recently, the BaBar collaboration reported their measurements for the ratios of 
the branching ratios of the semileptonic decay B — > D^r~u T and B — > D^*'l~Vi with 
i" = (e- AT) Q: 

g(i? -> Dr-P ) = 
V ; B{B^Dl-ui) 

where the isospin symmetry relations TZ(D°) = TZ(D + ) = TZ(D) and TZ(D*°) = TZ(D* + ) = 
1Z(D*) have been imposed, and the statistical and systematic uncertainties have been 
combined in quadrature. These BaBar results are surprisingly larger than the standard 
model (SM) predictions j^-S]: 

K{D) SM = 0.296 ± 0.016 [2], K(D*) SM = 0.252 ± 0.003 [2], 
K{D) SM = 0.316 ±0.014 [4], 

TZ(D) SM = 0.315 [5], K{D*) SM = 0.260 [5]. (2) 



The combined BaBar results disagree with the SM predictions by 3.2cx [4] or 3.4cr [1[. 
Since the report of BaBar measurements, this R(D^) anomaly has been studied in- 



tensively by many authors, for example, in Refs. [4 T2] - Some authors treat this 3.4er 
deviation as the first evidence for new physics (NP) in semileptonic B meson decays 
to t lepton 0-10], such as the NP contributions from the charged Higgs bosons in the 



Two-Higgs-doublet models 0] ■ Some other physicists, however, try to revisit the SM pre- 
diction for R(D^) by employing the relativistic quark model j5[, by using the form factors 
computed in unquenched lattice QCD or by maximally employing the experimental 
information on the relevant form factors from the data of B — > Dli>\ with I = (e~,/i~) 
[l3T |. but they all failed to explain the BaBar's R(D^) anomaly. 

In this paper, we will calculate the decay rates of the semileptonic B — > D^*'l~9i decays 
and the six R(X)-ratios: the four isospin- unconstrained ratios R(D°) ,R(D*°) ,R(D + ) and 
R(D* + ), as well as the two isospin-constrained ratios R(D) and R(D*) in the framework 
of the SM by employing the pQCD factorization approach (l4l |. We will compare the 
pQCD predictions with other theoretical predictions jH, in the framework of the SM, 
and the measured values of BaBar Collaboration jl[. 

In the pQCD approach, the lowest order Feynman diagrams for B — > D^l^ui decays 
are displayed in FigJTJ In the B meson rest frame, we define the B meson momentum Pi, 
the momentum P2 in the light-cone coordinates as 

Pr = ^|(l,l,0j, P 2 = ^V,rr,0 ± ), (3) 

where the ratio r = m D /m B or m D */m B , while = —r/~, 0±) / a/2 and €t = (0, 0, 1) 
denotes the longitudinal and transverse polarization of the D* meson respectively. The 

factors r/ ± = i] ± \/r\ 2 — 1 is defined in terms of the parameter 77 = i 1 + r 

The momenta of the spectator quarks in B and D^*> mesons are parameterized as k\ 
(0,Xi^,/ci_l) and k 2 = {x 2 rr] + , x 2 rr]- , k 2 ±). 
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FIG. 1. The lowest order Feynman diagrams for the semileptonic decays B 
pQCD approach, the winding curves are gluons. 



D^l'vi in the 



For the B meson wave function, we make use of the same one as being used in Refs. (l6k 
18| . We adopted the B-meson distribution amplitude widely used in the pQCD approach 



b) = N B x 2 (l — a;) 2 exp 



Mix 2 1, ,,„' 



(4) 



where Nb is the normalization factor, f B is decay constant. For parameter ug, one usually 
take ujb = 0.4 ± 0.04 GeV for B° and B mesons based on currently available data. 

For the pseudoscalar D meson and the vector D* meson, their wave function can be 
chosen as 19| 



$ D (j>, x) = -^=7 5 (^d + m D )(f) D (x) 



(5) 



® D *(p,x) = -J= \tLtyD*+m D .)<fe.(x)+t T tyD* + m D .)<fl ) .{x)\ . (6) 



For the distribution amplitudes of meson, we adopt the one as defined in Ref. 19 



2^6 



6x(l -x)[l + C DW (1 - 2x)\ exp 



cu 2 b 2 



(7) 



From the heavy quark limit, we here assume that fp* = /^* = fn*, 4>d* = 0d* = 4>d 



and set C D = C D . = 0.5, u = 0.1 GeV as Ref. 19 



The form factors Fo i+ (g 2 ) for B — >■ D transitions are defined in the same way as in 
Ref.[18|. In terms of fi(q 2 ) and fi(q 2 ), they can be written as 



F + (q 2 ) = -[f\(q 2 ) + f2(q< 



F (q 2 ) = -h(q 2 ) 



1 + 



+ \f2{q 2 ) 



1 - 



(8) 



where q = pi — P2 is the momentum of the lepton pair. Following the same procedure as in 
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Ref. 18], we calculate and obtain the B — >■ D form factors f\{q 2 ) and f 2 {q 2 ) as following, 
/i(<? 2 ) = 87Tm 2 B C F J dxidx 2 J bidb 1 b 2 db 2 <pB(xi, h)<f>D(x 2 , h) 

:{[2r (1 - rx 2 )) ■ h 1 (xi,x 2 ,b 1 ,b 2 ) ■ a s (ti) • exp [-S ab (tt)) 



x 



+ 



2r(2r c - r) + x x r -2 + 2r] + a/// 2 - 1 



2rj 



Vv 2 



+ 



v 7 ^ 



•/i 2 (^i, x 2 , h, b 2 ) ■ a s (t 2 ) ■ exp [-S ab (t 2 )] 



(9) 



f 2 (q 2 ) = Mm 2 B C F J dx x dx 2 J hdbib 2 db 2 (f) B (x 1 ,bi)(j) D (x 2 ,b 2 ) 
x \ [2 - 4x 2 r(l - rj)] ■ hi(x!, x 2 , b 1: b 2 ) ■ a s (t x ) ■ exp [-S ah {t x )] 



+ 



4r — 2r c — x\ + 



.1'! 



1 



h 2 (x 1 ,x 2 ,bi,b 2 ) ■ a s (t 2 ) ■ exp [-S^fo)] k (10) 



where CV = 4/3isa color factor, r c = m c /mB with m c is the mass of c-quark. The hard 
functions hi j2 (xi,bi) come form the Fourier transform and can be written as 

h 1 (x 1 ,x 2 ,b 1 ,b 2 ) = K (/3ibi) {6(b x - 6 2 )4(ai&2)^o(«i&i) 
+6(b 2 - bMa^K^h)} ■ S t (x 2 ), 

h 2 (xt, x 2 , h, b 2 ) = Ko^h) {9(bi - b 2 )I (a 2 b 2 )K (a 2 bi) 
+0{b 2 - bi)I {a 2 bi)K {a 2 b 2 )} • S t (x 2 ), 



'ID 



where K and Jo are modified Bessel functions, while the parameters «i = iriB^x 2 rr] + , 
a 2 = rriB\/xirr] + — r 2 + r 2 and 0x = f3 2 = mn\/x^x<?rr] + with r = m D {*) /tub- The 
threshold resummation factor St(xj) is adopted from [20]. The Sudakov factor S a b(t) = 



S B {t) + SM{t) can be written as |20j, \21 



s B (t) 



m B \ [ l djj, 



S M (t) = s(x 2 ^rri + M^+s(^l-x 2 )^r V + M^+2^ ^(a.^)), (12) 

with the quark anomalous dimension 7 g = —a s /ir. The explicit expressions of the func- 
tions s(Q,b) can be found for example in Appendix A of Ref. j2l[. The hard scales in 
Eqs.f l9lfT0|T2|) are chosen as the largest scale of the virtuality of the internal particles in 
the hard 6-quark decay diagram, 



ti = max{mBV 'x 2 rr] + , l/&i, I/&2}, h = ma,x{m B ^/x 1 rr] Jr - r 2 + r 2 , l/&i, 1/6 2 }.(13) 

With F + (q 2 ) and F n (q 2 \ the differential decay widths of the semileptonic decays £> — >• 
DWi can be written as 122] 



dT(b -> c/P,) c||y c 



dg 2 



1927T 3 



13 



1 



m 



2\ 2 



2g 2 

+ (m 2 + 2q 2 ) X(q 2 )\F + (q 2 )\ 2 }., 



3m 2 (m| -m 2 D ) 2 \F (q 



,2m2 



(14) 



where mi is the mass of the charged leptons, and A(g 2 ) = (m 2 B + m 2 D — q 2 ) 2 — Am 2 B m 2 D is 
the phase space factor. 

For the B — >■ D* transitions, the hadronic matrix elements of the vector and axial- 
vector currents are described by four QCD form factors V(q 2 ) and A 0t i^(q 2 ) [22]. By 
employing the pQCD approach, we calculate and find the form factors V(q 2 ) and Ao^^q 2 ) 
for B D* transitions: 

V (q 2 ) = Mm 2 B C F J dx x dx 2 J b 1 db 1 b 2 db 2 (f>B{xi, bi)4>D*( x 2, h) • (1 + r) 
x \ [1 - rx 2 ] ■ h 1 (x 1 ,x 2 , &!, b 2 ) ■ a s (*i) • exp [-^(ti)] 



+ 



r + 



h 2 (x 1 ,x 2 ,b 1 ,b 2 ) ■ a s {t 2 ) ■ exp [-S ab (t 2 )]j . 



(15) 



Ao(q 2 ) = 8Tcm 2 B C F J dx x dx 2 J b\db\b 2 db 2 (j)B{x\,bi)4' 1 D*{x 2 ,b 2 ) 

x j[l + r - rx 2 (2 + r - 2r?)] • hi^i, ar 2 , &i, &a) ■ a a (*i) ■ exp [-S ab {ti)\ 

= (l - 277(77 + V»7 2 - 1)) 



2 Xi 

r + r c + — + 



2^ - 1 2y/r] 2 - 
■h 2 (x!, x 2 , b h b 2 ) ■ a s (t 2 ) ■ exp [S ab (t 2 )] 



(16) 



Ma 2 ) 



87rm 2 B CF J dx\dx 2 J b 1 db 1 b 2 db 2 (p B (xi, bi)(f)J ) »(x 2 , b 2 ) ■ jq— 



x|2[l + r/ - 2rx 2 + rr)x 2 ] ■ h 1 (x 1 , x 2 , h, b 2 ) ■ a s {tx) ■ exp[-5 ab (ti)] 

+ [2r c + 2r/r - x x ] ■ h 2 (x 1 ,x 2 ,bi,b 2 ) ■ a s (t 2 ) ■ exp[-S ab (t 2 )} \, (17) 



2 _ (l + r) 2 (r/-r) 
^ } " 2r(77 2 - 1) 



Ax(g 2 ) — 8T[m 2 B Cp J dx\dx 2 J bidbib 2 db 2 (j) B {xi,bi) 



•0^(x 2 , b 2 ) ■ x { [(1 + 77XI - r) - rx 2 {l - 2r + 77(2 + r - 2r/))] 

•/ii(xi, x 2 , 61, 62) • a s (*i) • exp \-S ah {t{)\ 



+ 



X\ 

r + r c (r) — r) — r\r 2 + rxiff (77 + r) + x\ I r\r 



h 2 (x t ,x 2 , bi, b 2 ) ■ a s {t 2 ) exp [-S ab {t 2 )\ 



(If 



where r = mo* /tub, while Cp and r c is the same as in Eqs. (l9lfT0l) . And the hard function 
h(xi,bi), hard scales ti j2 and Sudakov factors S ab (ti) have been given in Eqs. (jlltil3p . 
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For B — > D*lv% decays, the differential decay widths can be written as 23 
dV L (B ->• D*W t ) G 2 F \V cb \ 2 



dq 2 



192vr 3 m| 



^m^X^q 2 ) 



2q 2 



mf + 2q 2 
Am 2 



(m B — m — q )(m B + m)Ai(q ) 



3m 2 X(q 2 )A 2 (q 2 
X(q^ 



m,B + m 



AM 



; (19) 



dT±(B D*Wi) 
dq 2 



G 2 F \V cb \ 2 
f927r 3 m| 



x 



f - 



m. 



[m] + 2q 2 



A 3/2( g 2 



V(q 2 



rxiB + rn 



(m B + m)A 1 (q 2 

7W) 



(20) 



where m = mo., and X(q ) = {m B + m|,* — g ) — Am B m D * is the phase space factor. 
The combined transverse and total differential decay widths are defined as 



dTn 



<IV.,. (IV. 



dq 2 dq 2 dq 2 



dV 

dq 2 



dV t dV T 
— - H 

dq 2 dq 2 



(21) 



In the numerical calculations we use the following input parameters (here masses and 

1.870, m D *o = 2.007, 



10, m 



decay constants in units of GeV) [18j, 
m D *+ = 2.010, tub = 5.28, m r = 1.777, m c 
r B ± = 1.641 ps, t b o = 1.519 ps, \Vcb\ = (41.21 
fD\/m D /m D *. 

By using the expressions as given in Eqs. 



mijo = 1.865, m£,+ 



1.35 ± 0.03, f B = 0.21 ± 0.02, f D = 0.223, 



0.5) 



x 10" 3 , A 



(/=4) 
MS 



0.287 and f D * 



and the definitions in Eq. 



we calculate the values of the form factors F 0:+ (q ), V(q ) and y4 0il 2 (g ) for given value 
of q 2 in the range of mf < q 2 < (itlb — m) 2 with m = vtld or m^*). Since the pQCD 
predictions for the considered form factors are reliable only for small values of q 2 (i.e. the 
large recoil regime), one has to make an extrapolation for them from the lower q 2 region 
to large q 2 region. In this work we make the pQCD calculations for the form factors in 
the lower q 2 
m 2 < q 2 < [tub 



region m" e < q" 

m) 2 by using the formula as given in Ref. [2i 



< m 2 , and then make the extrapolation to the range of 



F{q 2 



F(0) 



I -a- 



2 ■ 



(22) 



where F stands for the form factors F 0j+ ,V and ^4o,i,2, and a, b are the parameters to 
be determined by the fitting procedure. In Table QJ we list the pQCD predictions for all 
relevant form factors for B — )■ D and B — )• D* transitions at the scale q 2 = and q 2 = m 2 , 
respectively. The theoretical error of the form factors as shown in Table [I] is the total 
error: a combination of the two major theoretical errors from u>b = 0.40 ± 0.04 GeV and 
f B = 0.21 ± 0.02 GeV. 

After making numerical integrations, we find the pQCD predictions for the branching 
ratios of the eight B — > D^*H~Ui decays as listed in Table HTl where the individual the- 
oretical errors have been added in quadrature. In order to show the relative size of the 
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TABLE I. The pQCD predictions for the form factors of the B — > (D, D*) transitions at lower 
q 2 region for q 2 = and q 2 = m 2 . 





F (0) F+(0) FoK) *+K) 


B -> D 


0.561°;!? 0.56i°i? 0.69±°;il 0.75l°;!J 




C(0) 4,(0) Ai(0) A 2 (0) 


B -> D* 


U.Dd_ 012 U.OU_ 009 U.OZ_ 01Q U.O4_ 01Q 




V(m 2 T ) A (m 2 T ) Ai(m 2 T ) A 2 (m 2 T ) 


B -> D* 


0.8518:8 0.661°;!? 0.62^ 0.711^ 



errors, we here present the pQCD predictions for B(B — > D°'*°t u t ) with the errors: 

B(B~ -> D°r-P r ) = (117l^(a; fl )i°;i(/ s )l°^(y c6 )iHiK))%, (23) 
B(2T -> D*°r-P T ) = (1.82«;l(^)l°;l(/ B )l^(K b )lHoK))%, (24) 

where the four major theoretical errors come from the uncertainties of the input pa- 
rameters u B = 0.40 ± 0.04 GeV, f B = 0.21 ± 0.02 GeV, |V c6 | = (41.2l^) x 10" 3 and 
m c = 1.35 ± 0.03 GeV. From the numerical results in Table Ull where the errors of the 
pQCD predictions have been added in quadrature, one can see that since the theoret- 
ical errors are still large, around 35% in size, the pQCD predictions for the branching 
ratios agree with the measured values as given by BaBar Collaboration |l[ or quoted from 



PDG-2012 24 within one standard deviation. 



TABLE II. The pQCD predictions for the branching ratios of B — > D^*'t Ut and B -> D^l~ 
decays where l~ stands for e~ or j£~ . The world averages from PDG 2012 24|, and the BaBar 



measured values as given in Ref . jl , [25j, [26j] 



Channel 



P QCD(%) PDG(%)[24] BaBar(%)[l, 25, 26] 



B°- 


* D+t- 


v T 


1.07 


f0.44 
-0.34 


1.1 ±0.4 


1.01 ±0.22 


B°- 


> D+r 




2.51 


+-1.17 
-0.86 


2.18 ±0.12 


2.15 ±0.08 


B- - 




V T 


1.17 


f0.47 
-0.38 


0.77 ±0.25 


0.99 ±0.23 


B- - 




n 


2.7i; 


fl.26 
-0.93 


2.26 ±0.11 


2.34 ±0.14 


B°- 


-> D* + T 




1.67 


f0.48 
-0.46 


1.5 ±0.5 


1.74 ±0.23 


B°- 


* D*+r 


~n 


5.56 


fl.83 
-1.62 


4.95 ±0.11 


4.69 ± 0.34 


B- - 


-> D*°T 




1.82 


f0.52 
-0.50 


2.04 ±0.30 


1.71 ±0.21 


B- - 


-»■ D*°l- 


~n 


6.03 


fl.97 
-1.75 


5.70 ±0.19 


5.40 ± 0.22 



Since the theoretical errors of the pQCD predictions will be greatly canceled in the 
ratios of the corresponding branching ratios, we also define the six R(X)-ratios in the 
same way as Ref. [l| and compare our pQCD predictions with those reported by the 
BaBar measurements [lj]. For the two isospin-constrained ratios R{D) and R(D*), for 
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example, we find 



R{D) = B(B° D + l~u{) + B(B~ — >■ D°l~ui) ~ " 429 W^-o.orsK), (25) 

where the major theoretical errors come from the uncertainties of u>b = 0.40 ± 0.04 
GeV and m c = 1.35 ± 0.03 GeV. The theoretical errors of the pQCD predictions for 
the branching ratios induced by the variations of fs and \V c f,\ are canceled completely 
in the ratios of the branching ratios. It is easy to see that the theoretical errors of the 
pQCD predictions for R(X)-ratios are reduced significantly when compared with those for 
branching ratios, from around 35% to ~ 5% only. 



TABLE III. The pQCD predictions for the six R(X)-ratios, other SM predictions @, 0, S] , and 
the measured values from BaBar Collaboration 111]. 



Ratio 


pQCD 


SM [2j SMf4] SM[5] 


BaBar [1] 


R(D°) 
R(D+) 
R(D*°) 
R(D* + ) 


n 4 o 2 +0.021 
U.UZ_q Q24 

U.4ZD_ 22 

n 302+ - 011 

u.ouz_ 016 
u.ouu_ 015 




0.429 ± 0.097 
0.469 ± 0.099 
0.322 ± 0.039 
0.355 ± 0.044 


K(D) 
K(D*) 


0.429 ± 0.023 

S01+ 012 
u.oui_ 015 


0.296 ±0.016 0.316 ±0.014 0.315 
0.252 ± 0.003 - - 0.260 


0.440 ± 0.072 
0.332 ± 0.030 



In Table IIII1 we list the pQCD predictions for all six R(X)-ratios, the SM predictions 
as given in Refs. 0, 0, [EJ respectively, and the measured values as reported by BaBar 
Collaboration [l[. From the numerical results one can see that 

(i) Due to the strong cancelation of the theoretical errors in the ratios of the corre- 
sponding branching ratios, the error of the pQCD predictions for both R(D) and 
R(D*) is ~ 5% only, similar in size with other two SM predictions 0, 0]. 

(ii) In Refs. 0, [EJ , the authors calculated the ratios by using the hadronic form factors 
computed in unquenched lattice QCD[4| or by employing the relativistic quark 
modelj^j, but as shown by the numerical results in Table IHIt they all failed to 
explain the BaBar's anomaly of R(D(*)). 

(hi) For the four isospin-unconstrained ratios R(D°) , R(D + ) , R(D*°) and R(D* + ), the 
pQCD predictions agree well with the BaBar measured values within one standard 
deviation. 

(iv) For the two isospin-constrained ratios R{D) and R(D*), the pQCD predictions are 
R(D) = 0.429 ± 0.023 and R(D*) = 0.301±g;gi§, the BaBar's anomaly of R(D^) 
are therefore explained successfully in the framework of the SM by employing the 
pQCD factorization approach. 
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In summary, we studied the semileptonic decays B — > D^*'l~ui in the framework of the 
SM by employing the pQCD factorization approach. We firstly evaluated the relevant form 
factors F 0i+ (g 2 ), V(q 2 ) and t4 0i i i 2(<? 2 ), and then calculated the branching ratios B(B — > 
D^*>l~u{) and the six R(X) ratios. From the numerical calculations and phenomenological 
analysis we found that 

(i) The pQCD predictions for B(B ->• D^l'vi) agree well with the measured values 
within one standard deviation. 

(ii) For the isospin-constrained ratios R{D) and R(D*), the pQCD predictions are 
R(D) = 0.429 ±0.023 and R(D*) = 0.301^^, the BaBar's discrepancy of 2.0a for 
R{D) and 2.7a for R(D*) are reduced significantly to 0.2a and 1.0a, respectively. 
We provide successfully the first SM interpretation for the BaBar's R(D^) anomaly 
by employing the pQCD approach. 
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